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THE  GAS  FLOW  IN  GAS-OPERATED  WEAPONS 


ABSTRACT 

In  gas-operated  veapons,  the  time-varying  pressure  in  the  barrel  is 
fed  through  a  duct  into  a  cylinder.  The  piston  in  the  cylinder  is 
displaced  by  the  pressure  and  operates  on  a  mechanism  which  extracts 
the  spent  cartridge  and  completes  the  next  loading  cycle. 

f 

The  theory  presented  here  predicts  the  pressure  history  in  the 
cylinder  and  the  motion  of  the  piston  for  a  given. pres sure  and 
temperature  history  in  the  barrel. 
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I .  INTRODUCTION 


In  gas  operated  weapon#  or  in  control  system#  a  high  pressure,  high 
temperature  supply  of  gas  is  used  to  operate  certain  mechanism#.  The 
crucial  problem  presents  itself  in  the  following  form:  At  a  given  point, 
A,  in  the  system,  pressure  and  tempera ture  are  prescribed  as  a  function 
of  time  due  to  the  inherent  operation  characteristics  of  the  system. 

This  pressure  is  communicated  through  a  duct  to  another  point,  B,  in  the 
system  where  the  gas  is  used  to  operate  a  mechanism;  for  example,  extract 
a  spent  cartridge  and  impart  enough  momentum  to  the  bolt  so  that  the  next 
round  may  be  positioned  in  the  chamber.  It  is  desired  to  predict  the 
pressure  and  other  variables  at  point  B  as  a  function  of  time.  The  mag¬ 
nitude  and  time  history  of  this  pressure  can  be  crucial  to  the  whole 
operation  cycle.  In  the  example  cited  above,  the  pressure  pulse  rust 
not  cause  extraction  at  too  early  a  stage,  when  the  cartridge  is  still 
pressed  to  the  wall  of  the  chamber  by  the  pressure  in  the  chamber;  on  the 
other  hand  the  pressure  history  must  be  such  that  sufficient  momentum  is 
transferred  to  the  bolt  to  operate  the  weapon.  The  pressure  history  is 
coupled  to  the  rate  at  which  work  is  extracted  at  B  from  the  gas  and  also 
to  the  nature  of  the  flow  of  gas  between  A  and  B.  In  this  report  an 
analysis  is  developed  which"  describes  the  time  history  of  pressure  and 
temperature  at  point  B  for  prescribed  condition  at  A  in  the  practically 
important  case  wl  re  the  characteristic  time  of  the  pressure  pulse  at  B 
is  large  compared  to  the  time  an  acoustic  signal  needs  to  travel  from 
A  to  B. 

The  analysis  is  developed  on  the  specific  example  of  the  M-l6  rifle 
where  time  varying  pressure  is  fed  from  a  port  in  the  barrel  (point  A)  to 
the  bolt  cavity  (point  B)  where  the  pressure  is  used  to  unlock  the  bolt, 
extract  the  cartridge  and  supply  sufficient  momentum  to  the  bolt  carrier 
so  that  the  next  loading  cycle  can  be  completed.  The  analysis  applies, 
however,  equally  well  to  other  gas  operated  weapons  and  should  find 
application  in  rocket  control  systems  where  high  pressure  gas  is  used  to 
operate  jet  or  surface  controls. 


15 


Preceding  Page  Blank 


II.  METHOD  OF  APFKOACH 


For  purposes  of  orientation.  Figure  1  may  be  consulted .  This  shows 
a  design  drawing  of  the  gas  system  of  the  M-l6  rifle.  As  mentioned 
before,  hot,  high  pressure  powder  gases  are  extracted  from  the  barrel  and 
fed  into  the  bolt  chamber  by  a  duct.  The  pressure  variation  with  time  at 
the  port  is  determined  by  the  gas  flow  in  the  barrel.  The  mass  flow 
through  the  port  is  small  compared  to  the  mass  flow  in  the  barrel  itself, 
and  the  effect  of  mass  removal  is  negligible  on  the  pressure  history  at 
the  port.  Thus,  the  pressure  and  temperature  may  be  considered  prescribed 
and  known,  say  from  an  interior  ballistics  analysis  of  the  weapon.  The 
hot  gases  expand  in  the  bolt  cavity,  thereby  accelerating  the  bolt  carrier, 
C.  Unlocking  of  the  holt,  B,  is  accomplished  after  the  carrier  has 
traveled  a  certain  distance,  at  which  distance  vent  holes  are  laid  free 
and  the  gas  in  the  bolt  carrier  expands  to  atmosphere.  At  this  distance 
the  holt  carrier  engages  the  bolt  and  cartridge,  and  extraction  is 
accomplished  by  transferral  of  energy  and  momentum  of  the  carrier  to  the 
bolt  and  cartridge.  The  residual  momentum  and  energy  of  the  bolt  carrier 
is  sufficient  to  complete  the  cycle  of  cocking  the  weapon  and  positioning 
a  new  cartridge  in  the  chamber.  It  is  the  object  here  to  predict  for  a 
given  pressure  and  temperature  at  the  port  the  motion  of  the  bolt  carrier 
referenced  in  time  to  the  rise  of  pressure  at  the  port.  The  motion  of 
the  bolt  carrier  is  governed  by  the  pressure  and  temperature  in  the  bolt 
carrier  cavity. 

It  is  natural  to  divide  the  analysis  into  twc  parts.  The  first  is 
the  determination  of  the  "space  averaged"  pressure  and  temperature  in 
the  bolt  carrier  cavity  for  given  flow  conditions  at  the  entrance  of  the 
cavity.  In  general,  these  depend  on  the  conditions  in  the  cavity.  It 
will  be  shown  in  Section  III. A  that  conservation  of  energy  and  mass 
together  with  the  equation  of  motion  of  the  carrier,  which  here  takes  the 
place  of  conservation  of  momentum,  uniquely  determine  the  conditions  in 
the  cavity  and  the  motion  of  the  carrier.  The  second  part  considers  the 
flow  from  the  port  to  cavity.  The  unsteady  flow  in  the  pipe  is  quite 
complex  and  is  indeed  not  even  amenable  to  numerical  computations  without 


1 6 


C  (8oit  Carrier) 
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B(8olt  Carrier  Cavity) 


considerable  simplifications.  Such  a  simplification  is  offered  by  consid¬ 
ering  the  flow  to  be  one- dimensional;  i.e.,  flow  variables  are  assumed 
constant  over  the  cross  section  of  the  duct,  so  that  the  flow  variables 
depend  only  cn  x,  the  distance  along  the  axis,  and  on  time*  Implicitly, 
this  assumes  that  conversion  of  frictional  work  to  heat  and  conduction 
of  heat  occur  instantaneously  and  uniformly  over  the  entire  cross  section. 

It  is  usual  to  introduce  friction  (and  heat  transfer)  even  in  time 
dependent  compressible  flow  according  to  the  practice  in  hydraulics  by 
means  of  a  friction  coefficient  f  defined  by  F  =  2  f  u,2/d  where  F  is  the 
friction  force  per  unit  mass.  The  friction  coefficient  f,  however,  is 
not  known  in  unsteady  turbulent  flow  and  its  prediction  would  indeed 
require  the  solution  of  the  full  problem  since  the  coefficient  of  friction 
will  depend  on  the  entire  flow  history.  Customarily,  the  friction  co¬ 
efficient  for  steady  flow  is,  therefore,  used  and  the  few  experiments 
done  indicate  this  to  be  a  good  approximation1  as  long  as  no  separation 
occurs.  However,  it  is  well  to  keep  in  mind  that  this  assumption  is 
satisfactory  only  for  nearly  steady  flow. 

It  will  be  shown  in  Section  III.C  that  the  flow  here  can  be  reduced 
to  a  quasisteady  flow  except  for  the  very  early  phases  when  the  flow  is 
being  established.  For  the  quasisteady  regime,  the  frictional  effects 
may  be  justifiably  included  by  using  the  steady  friction  coefficient. 

For  this  phase  of  the  flow,  pressure  and  momentum  losses  due  to  sudden 
enlargements  or  constrictions  are  included  in  the  analysis  within  the 
framework  of  one- dimensional  quasisteady  flow.  Losses  due  to  the  various 
bends  in  the  duct  can,  for  large  -t/D  of  the  duct,  be  included  by  loss 
coefficients.  These  losses  will  be  treated  as  a  distributed  resistance 
and  accounted  for  in  an  overall  friction  coefficient. 

For  the  very  early  stages,  the  flow  is  truly  unsteady  and  it  is  not 
possible  to  include  the  frictional  effects  using  a  steady  friction 
coefficient,  nor  is  it  possible  to  consider  the  bends  in  the  tube  as 
distributed  resistances.  This  early  stage  involves  a  strong  shock, 
which  actually  starts  the  flow;  this  shock  will  be  reflected  and  dif¬ 
fracted  and  also  influenced  by  friction  and  heat  transfer  in  a  complex 


I 

:i 

way  not  accessible  to  analysis  at  this  time.  The  first  stage  will,  ; 

therefore,  he  treated  entirely  without  friction  and  discarding  all 

possible  reflections,  the  only  justification  being  that  this  first  stage 

is  of  very  short  duration  compared  to  the  total  flow  time.  As  will  be 

shown  in  Section  III.D  ,  the  first  stage  essentially  introduces  a  delay 

time  between  the  pressure  rise  at  the  port  and  the  pressure  rise  in  the 

cavity.  The  net  effect  of  neglecting  friction  and  shock  reflections  in 

the  first  stage  will  be  to  predict  a  delay  time  somewhat  too  small. 

However,  this  delay  time  has  been  observed  experimentally  in  numerous 

p 

tests  and  a  correction  to  the  theoretical  value  is  possible,  thus 
accounting  for  the  observable  effect  of  friction  and  shock  reflections 
in  the  first  stage.  The  gas  wi]l  be  considered  a  calorically  and 
thermally  perfect  gas ;  thus ,  the  ratio  of  specific  heats  and  the  gas 
constant  will  be  considered  constant  and  given,  say  from  the  interior 
ballistic  analysis  of  the  weapon.  Because  of  the  high  pressure  and 
density  the  assumption  of  equilibrium  for  the  gas  phase  of  the  propel¬ 
lant  can  be  safely  made;  however,  there  is  nonequilibrium  between  the 
solid  and  the  gas  phase  of  the  propellant  in  the  barrel  itself  and  also 
in  the  duct.  Within  certain  temperature  and  density  ranges  it  is  always 
possible  to  approximate  the  real  gas  behavior  by  an  ideal  gas  behavior 
as  long  as  the  gas  phase  is  in  thermodynamic  equilibrium.  This  approxi¬ 
mation  is  commonly  made  in  the  standard  interior  ballistic  treatment, 
and  performance  predictions  based  on  this  assumption  are  apparently  in 
good  agreement  with  experiment.  Therefore,  assumption  of  equilibrium  and 
indeed  perfect  gas  behavior  is  expected  to  be  adequate  also,  since  the 
primary  purpose  here  is  to  clear  up  the  fluid  dynamic  aspects.  Correc¬ 
tion  factors  accounting  for  real  gas  effects  may  be  applied  to  the  re¬ 
sults  found  on  the  basis  of  perfect  gas  behavior  as  shown,  for  example, 
in  reference  3.  It  should  be  mentioned  that  a  rigorous  nonequilibrium 
treatment  of  the  gas  flow  is  not  possible  because  the  kinetics  of 
powder  gas  reactions  are  not  sufficiently  understood. 

According  to  the  assumption  made  above,  the  analysis  will  be  based 
on  the  simplified  geometry  shown  in  Figure  2.  The  area  restriction  at 
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Figure  2.  Idealized  Geometry  for  M-l6  Rifle 


the  port  will  be  treated  as  a  supersonic  nozzle  of  zero  length  when  the 
flow  into  the  duct  is  supersonic, as  is  indicated  on  Figure  2.*  The  duct 
is  considered  a  straight  tube  of  constant  inner  cross  section,  with  all 
losses  resulting  from  bends  and  sudden  area  changes  taken  into  account 
as  discussed  above. 

Heat  transfer  from  the  gas  to  the  wall  is  assumed  to  leave  the  wall 
temperature  unchanged;  this  assumption  allows  uncoupling  of  the  fluid 
dynamic  problem  from  the  heat  conduction  problem  in  the  duct,  &nd  is  a 
good  one  for  the  small  Fourier-numbers  of  interest  here. 

III.  ANALYSIS 


A.  The  Flow  in  the  Cavity 

The  average  properties  in  the  cavity  are  determined  by  the  conserva¬ 
tion  equations  of  mass  and  energy  and  by  the  equation  of  motion  of  the 
bolt  carrier. 

Let  in  =  Ag  pe  ug  (l) 

be  the  mass  flow  into  the  cavity,  where  A  is  an  area-x-^p  the  gas  density 
and  the  ii  the  gas  velocity,  and  where  the  subscript  designates  the  loca¬ 
tion  of  these  quantities  in  accordance  with  Figure  2.  Conservation  of 
mass  gives 


1  *  (Vcl  ♦  \  s„)  3!*  *  a,  p0  vB  ♦  K>  k.v  (v  pc  Pc)^  (2) 

The  first  term  on  the  right  side  is  due  to  the  change  of  density  in  time 
in  the  cavity;  the  second  term  is  due  to  the  change  in  cavity  size;  and 
the  third  term  is  the  term  due  to  venting. 

Accordingly,  =  0  for  xB  <  xgy  -rv,  and  Kj  =[2/(Y+l)lY+l)/2(Y'l) 
for  x.  >  Xgv  -  ry,  where  Xgy  and  rv  are  the  center  location  and  radius, 
respectively,  of  the  circular  vent. 

*The  case  where  the  port  is  treated  as  an  area  discontinuity  is  dis¬ 
cussed  in  the  Addendum. 

^Dimensional  quantities  are  indicated  by  bar. 


21 


Here,  it  has  been  assumed  that  is  always  much  larger  than  the 
ambient  pressure,  so  that  venting  occurs  with  sonic  velocity.  Mass 
losses  due  to  leakage  may  be  included  in  this  equation  by  an  obvious 
extension. 


Heat  transfer  to  the  wall  of  the  cavity  may  be  neglected,  since  the 
average  gas  velocity  in  the  cavity  is  small  so  uhat  convection  is 
negligible. 

Conservation  of  energy  may  then  be  written  as 

5  \  tot  =  [»c  So  V  2  5=  V  5o  Jl  h  a5]  + 

*1  ?o  5v  +  K2  Vo  (T  h'>c)1/2  *  <3) 

^vYPc  <Y5C/^)1/2 

The  terms  in  the  square  bracket  are  the  thermal  energy  of  the  gas  in  the 
cavity,  the  kinetic  energy  of  gas  in  the  cavity  and  the  work  done  by  the 
gas  in  the  cavity.  The  three  other  terms  are  the  energy  losses  due  to 
venting.  For  the  constant  Kg  we  have 

Kg  =  0  for  xB  <  xBv  -  ?v;  Kg  =  [2/(Y+l) ] (3y"1)/2(y_1)  for  ^  >  5^-  ?y 


For  the  circular  vent 

=  <V  V>  [?v2  -  <*B  -  *Zv)Z]112  *  \  [2  *  -  V)/?v)! 


for 


x„  -  r  <  <  x^.  +  r 

Bv  v  —  B  —  Bv  v 


T  -  2 

A  =  tt  r 
v  v 


for 


x  +  r  <  x_. 
3v  v  —  B 


The  volume  V  of  the  cavity  is  given  by: 

V  =  V  .  +  A  x_ 
c  ci  c  B 
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The  term  m  h  ,  .  is,  of  course,  the  flux  of  enthalpy  into  the  cavity, 

S  tOu 

where 

h.,=p/p+e+u/2 
e  tot  *e  e  e  e 

The  equation  of  motion  for  the  bolt  carrier  is  given  by: 


d  v  -1  -1 

-rmhlc"B  +$mb 

at 


(k) 


d^/dt  =  vB 

$  signifies  external  forces  on  the  bolt  carrier,  such  as  friction  on  the 
bolt.  At  the  end  of  the  bolt  carrier  travel  v^  and  dv^/dt  are  set  equal 
to  zero. 

The  conservation  equations  and  the  equation  of  motion  of  the  bolt 
are  three  equations  for  the  five  unknowns,  SL,  p  ,  p  »  6  and  e  . 

Q  v»  c  C 

The  problem  is  rendered  complete  by  the  addition  of  a  thermal  and 
a  caloric  equation  of  state,  thus: 


e  =  cv9 


and  p  =  R  p  0 

r  gas 


(5) 


B.  The  Flow  in  the  Duct 

The  flow  in  the  duct  is  a  mixed  initial  and  boundary  value  problem. 
Using  the  usual  notation  for  pressure,  velocity  and  sound  velocity,  p, 
u,  a,  respectively,  the  initial  conditions  are 

t  =  0;  x  >  0  :  p  (x,0)  =  p.  ; 


a  (x,0)  =  a^  ;  u  (x,0)  =  0 


x  <  0  :  p  (x,0)  =  p  (0);  a  (x,0)  =  a  (0);  u  (x,0)  =  0 
8  8 


(6) 


For  ^  >  Sp  the  boundary  x  =  0  is  spacelike  and  three  input  data  at  x  =  0 
determine  the  solution  for  nonisentropic  flow.  Thus, 

t  >  0,  x  =  0  :  p  (0,t)  =  p  (t);  a  (0,t)  =  a  (t)  ;  u  (0,t)  =  u  (t)  (7) 

r  x  r 


For  Up  <  the  line  x  =  0  is  timelike  corresponding  to  the  fact  that  only 

two  characteristics  enter  the  region  of  interest.  Thus ,two  data  along 
x  =  0  are  required: 

t  >  0,  x  =  0:  5(0, t)  =  PT3(t) ,  u(0,t)  =  UB(t)  (8) 

By  virtue  of  the  fact  that  the  nozzle  length  is  very  small,  the 
conditions  at  the  exit  of  the  nozzle  (subscript  p)  follow  the  conditions  in 
the  barrel  (subscript  g)  instantaneously. 

At  x  =  £-no  boundary  condition  can  be  prescribed  for  ug  >  ag 
corresponding  to  the  fact  that  no  characteristic  reaches  the  space  x  <  1 
for  dt  >  0.  For  ue  <  ag  the  characteristic  given  by  dx/dt  =  u  -  a 

reaches  into  the  region  x  <  1  and  one  condition  at  x  =  I  must  be  prescribed. 
Thus, 

t  >  0,  x  =  £:  p(£,t)  =  PQ(t)  (8a) 

Subject  to  the  above  initial  and  boundary  values  the  following  set  of 
equations  describes  the  one-dimensional  flow  in  the  duct. 


Continuity: 


Momentum: 


ii  +  ?J.P-  uj  =  o 

3E  3x 


+  i  |-2-  +  f  =  o 
D  l  p  3  x 


Energy : 


1  3p  ,  -  D  /U  ,  r\ 

- +  q  =  —  ( - +  h), 

p  3t  tit  2 


-  =  ,  -  De  p  Dp 

uF  +  q  = - -r- 

Dt  p  Dt 


(11a) 


Instead  of  the  energy  equation  the  entropy  equation  may  be  used: 

0  —  =  u  F  +  q  (.] 

Dt 

Here,  D/Dt  =  a/at  +  u  3/cx  is  the  material  derivative. 
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These  equations  are  complemented  by  the  equation  of  state  (5).  The 
hydrodynamic  friction  force  F  is  introduced  as 


F  = 


D  ir  t  dx 
w 

,«2 


I)  T  _o 

v  2  fu 


p  it  {u  Ik)  dx  p  D 


(13) 


D 


1 _ 2 

with  t  ®rf  p  u  ,  the  wall  shearing  stress. 

W  c. 


The  amount  of  heat  transferred  per  unit  time  per  unit  mass  is: 


q  =  - 


5(6  -6  )  D  77  dx  -1*5  (0  -  6  ) 
jS _ _ _  =  . g .  y.. 


(D  ^/k)  tt  p  dx 


(Hi) 


P  D 


Reynolds  analogy  for  turbulent  flow  gives: 

f  c  p  u 


J  = 


(15) 


This  relation  is  not  strictly  fulfilled  for  turbulent  flow  in  a  pipe 
because  of  the  effect  of  the  laminar  sublayer.  Here,  it  is  assumed  that 
this  relation  holds  strictly;  this  amounts  to  setting  the  recovery  factor 
equal  to  unity. 

The  above  system  of  equations  is  hyperbolic  and  has  the  character¬ 
istic  directions 

dx/dt  =  u  +  a  and  dx/dt  =  u  , 

the  latter  being  the  differential  equation  for  the  particle  path  and  the 
former  the  differential  equations  for  the  other  two  characteristic 
directions.  The  compatibility  conditions  may  be  found  by  standard 
techniques  and  are  given  in  many  textbooks,  see  for  example  reference  U. 
Thus  the  flow  in  the  duct  may  be  computed  by  the  method  of  character¬ 
istics.  However,  such  a  computation  for  the  present  case  is  very  time- 
consuming  at  best.  Since  the  acoustic  transit  time  I/a  is  very  short 
compared  to  the  total  duration  t  of  the  flow,  the  computation  has  to  be 
extended  to  large  times  and  this  will  lead  to  an  appreciable  error  unless 
the  grid  size  is  very  small.  It  should  be  pointed  out  that  character¬ 
istic  computation  for  nonisentropic  flow  is  considerably  raoi'e  complex 
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since  the  particle  path  has  to  be  computed  together  with  the  two  other 
characteristics.  A  further  complication  arises  from  the  fact  that  the 
boundary  conditions  are  such  that  the  solution  will  not  be  continuous 
throughout,  but  shocks  will  appear  imbedded  in  the  flow  field.  The 
method  of  characteristics  is  not  capable  of  handling  the  appearance  of 
shocks  without  appreciable  complications. 


The  analysis  proposed  here  specifically  exploits  the  fact  that  the 
acoustic  transit  time  is  small  acompared  to  the  operation  time  x.  This 
is  the  situation  in  gas  operated  systems  discussed  so  far;  and  the 
analysis  is  applicable  for  situations  where  the  numerical  computation  by 
the  method  of  characteristics  becomes  impractical. 


With  the  i  the  characteristic  length  dimension,  x  the  character¬ 
istic  time  and  a  (0)  the  characteristic  velocity  the  following  non- 
S 

dimensional  variables  are  introduced. * 


1)  t 

2)  x 

3)  u 

1+)  a 

5)  P 

6)  6 
7)  P 

3)  h 


tx 

xZ 

a  (0)u 

g 


a  (0)a 


P  (0)a  <0)p 

O  D 


a  2  (0)  e/R 


gas 


Pg  (0) p 


a  2(0)h 


9)  F  =  a  2(0)F/£ 
& 

10)  i  =  ag3(0)q/2 


(16) 


*The  sound  velocity  at  t  =  0  may  indeed  serve  as  reference  velocity 
since  the  sound  velocity  ag  (t)  is  not  a  strong  function  of  time 

a  ~  (p  )  W-1)/2Y  with  (y-1)/2y  4  1,  where  y  =-  1.25. 
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The  differential  equations  then  become 


.e  kSSL  +  eL 

p  at  at 


2fi.+ 

8x 

jlgal  =  o 

3x 

(17) 

+ 

at 

3u  .  1  3p  .  „ 
u  ~ —  +  -  -rs~  +  F  =  0 

3x  p  ax 

(18) 

2  a  2 

(h  + 

r )  - u  Ir (h  +  r>  - 4  =  0 

(19) 

where 


£HTV> 


is  the  ratio  of  acoustic  transit  time  to  the  flow  duration  time.  We 
note  here  that  the  flow  duration  time  7  is  of  the  same  order  as  the 


average  characteristic  time  p  /(9p  /at),  impressed  on  the  flow  in  the  duct 

&  D 

by  the  time  variation  of  the  pressure  in  the  barrel.  Thus,  we  have 
^  '  0  (pg/[a.pg/3t])  with  t  »£/ag(0) 

C.  The  Outer  Problem 

Since  e  is  a  small  number,  we  assume  asymptotic  expansions  for  the 
dependent  variables  of  the  form 

p  (x,t ;e)  =  p1  (x,t)  +  $2  (e)  pg  (x,t)  +  .  .  . 

P  (x,t;e)  =  p1  (x,t)  +  62  (e)  p2  (x,t)  +  .  .  . 

h  (x,t;e)  =  h1  (x,t)  +  $2  (e)  hg  (x,t)  +  .  .  .  (20) 

u  (x,t;e)  =  ux  ( x ,t )  +  $2  (e)  u2  (x,t)  +  .  .  . 

6n+l(e) 

here  6 — ID  *  °  in  the  limit  e  ^  0  811(1  <5-,  (O  =  1  since  the  leading 
n 

term  for  all  expansions  is  of  order  unity. 
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Substituting  these  expansions  into  the  full  actuations  ( IT )  >  (l8), 
and  (19)  gives  in  the  limit  e  -*■  0  the  following  set  for  the  first  order 
solution 


3  (piui) 
3x 


=  0 


3x 


3x 


+  F  =  0 


u. 


u. 


)  -  q  =  0 


(21) 


These  equations  are  the  equations  describing  the  quasisteady  flow; 
i.e.j  time  appears  in  these  equations  only  as  a  parameter.  Consequently, 
the  initial  condition  specified  by  (6)  cannot  be  fulfilled  by  the 
solution  to  the  above  equations  (2l).  This  is  a  result  of  loss  of  the 
terms  involving  time  derivative  by  the  expansion  process.  The  solution 
to  (21)  with  the  boundary  condition  (8)  and  (8a)  cannot  be  uniformly 
valid  in  t  because  the  initial  conditions  are  violated.  In  order  to 
obtain  a  uniformly  valid  solution  another  expansion  for  small  t  will  be 
needed.  This  expansion  will  be  given  in  Section  m.D. 

In  the  terminology  of  the  method  of  matched  asymptotic  expansions 
the  solution  to  (21)  is  called  the  outer  solution,  while  the  expansion 
for  small  t  to  be  developed  in  Section  III.  D  is  called  the  inner 
solution.  It  will  be  shown  later  that  both  solutions  can  be  matched 
and  a  composite  solution  can  be  found  which  is  then  a  single  expansion 
uniformly  valid  in  t. 

Customarily,  the  variables  pertaining  to  the  outer  solution  are 
designated  by  small  letters  and  those  pertaining  to  the  inner  solution 
by  capital,  letters.  The  subscript  "l",  designating  the  first  approxi¬ 
mation  in  the  outer  problem,  will  now  be  dropped. 
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The  outer  problem  is  that  of  a  viscous  heat  conducting  flow  in  a  pipe 
of  constant  diameter.  This  is  a  problem  of  eminent  practical  importance 
and  has  been  treated  as  early  as  1875  (reference  5).  Extensive 
discussions  may  be  found  in  references  3*  6,  and  7*  Here,  the  treatment 
given  in  reference  6  will  be  followed,  and  only  the  main  results  quoted. 
For  more  detailed  analysis  the  reader  is  referred  to  reference  6. 

Using  the  notation 


c  9,  .  =  h,  =  h  +  — 
p  tot  tot  2 


the  energy  equation  may  be  written  as 


-  q/c  =  q  R  /c 
p  gas  p 


q  =  -2f  -  (9,  .  -  6  )  u 

fl  D  tot  w 
gas 


we  have 


2ft 


(6tot  -  V 


Since  6^  is  assumed  constant  in  x  and  t,  this  equation  may  be  integrated 
with  the  boundary  condition  9+o^.  =  0^  at  x  =  0: 


0,  .  =  0  +  (0  -9  )  e 
tot  w  gw 


-2fIx/D 


Using  ( 21)  and  the  equation  of  state  (5)  one  obtains  with  M  =  u/a  the 
following  ordinary  differential  equation  (see  reference  6): 


(1  +  Y  M2)  (1  +  ^  M2)  d  etot  Y  M?(l  +  m2) 

M2  (1  -  M2)  °tot  ( 1  -  M2 ) 


U  f  -  dx  (2U) 
D 
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Together  with  equation  (22) ,  this  determines  the  Mach  number  distribution 
in  the  duct.  This  equation  cannot  be  integrated  in  closed  form  and  a 
numerical  solution  must  be  found.  Once  the  Mach  number,  M,  and  the  total 
temperature,  0^  &re  determined  the  other  dependent  variables  can  be 
found  quite  easily.  The  pertinent  relations  are  provided  later  on. 


First,  the  special  case  0^  ^  >>  0^  mey  be  considered.  In  this  case 
separation  of  variables  is  possible  if  0  is  neglected  in  comparison  with 
9tot*  ^ne  th8n  obtains  from  (2k)  and  (22)  the  differential  equation: 


a?  it  i  (i  +  W1 

2  D  1  -  M2 


This  equation  is  readily  integrated  in  closed  form,  though  the  solution 
is  omitted  here. 


From  the  differential  equation  one  notes  the  following  condition  for 
the  sign  of  dM2  for  supersonic  flow: 

M2  >  1  i  dM2  <  0 

This  means  the  Mach  number  always  decreases  if  the  initial  Mach  number 

at  x  =  0  is  supersonic.  Mach  number  one  is  the  limiting  Mach  number 

that  can  be  reached  on  the  supersonic  branch  and  this  Mach  number  is 

reached  for  a  distance  x  =  I  which  is  determined  by  the  friction 

coefficient  f,  but  which  is  quite  small.  If  the  duct  length  is  larger 

than  1  ,  then  a  shock  occurs  somewhere  in  the  duct  reducing  the  Mach 

max’ 

number  to  subsonic  values  and  the  subsequent  flow  stays  subsonic.  For 
the  subsonic  flow: 

-  <  M2  <  1 ,  dM2  >  0 
Y 

This  means  the  Mach  number  increases  in  the  direction  of  the  flow.  Hov~. 
ever,  M  =  1  cannot  be  surpassed  and  is  reached  only  at  the  exit  if  pg  > 
Pc>  since  dM2/dx  -*  «  as  M  *  1  so  that  the  inlet  conditions  are  adjusted 
to  allow  M  =  1  to  occur  only  at  the ’exit.  For  p  =  p  the  flow  is 

C  6 

subsonic  at  the  exit  also.  The  lower  limit  at  1/y  is  a  consequence  of 

2 

the  special  assumption  of  0w  =  0.  For  M  =  1/y  at  x  =  0,  the  Mach 


30 


number  stays  constant  in  the  duct  since  dM^/dx  =  0.  Below  this  value 

dM2  <  0,  i.e  . :  the  Mach  number  decreases.  As  equation  (23)  shows,  the 

condition  6  <<  9,  .  is  not  a  valid  assumption  for  the  specific  case  of 
w  tot 

the  M-16  rifle,  but  could  be  a  good  assumption  for  a  different  application. 

For  the  general  case  0  is  finite,  then  tne  Mach  number  in  the  subsonic 

branch  first  decreases  and  then  increases  reaching  the  condition  dM2 /dx  = 

0  only  for  one  value  of  x,  which,  of  course,  gives  the  minimum  value  of  M. 

M  =  1  is  here  reached  only  at  the  exit  for  p  >  p^;  for  p  -  p  only 

subsonic  flow  can  occur  at  the  exit.  The  supersonic  branch  has  the  same 

qualitative  behavior  as  above.  The  phenomenon  discussed  above  is  termed 

"choking"  because  it  gives  a  limitation  on  the  maximum  mass  flow.  The 

conditions  for  choking  in  the  general  case  are  difficult  to  state  but  a 

7 

lucid  physical  interpretation  is  given  by  Prandtl  for  the  case  of  zero 
heat  transfer. 

Figure  3  shows  the  Mach  number  distribution  along  x  from  numerical 

computations  for  the  general  case  of  0  finite.  Shown  also  is  the 

w 

location  of  the  shock.  For  a  given  Mach  number  distribution  in  the  super¬ 
sonic  branch,  say  M^  =  f(x),  and  a  given  Mach  number  distribution  for  the 
subsonic  branch,  say  M^  =  g(x) ,  the  location  of  the  shock  xg  is  uniquely 
determined  by  the  condition: 


For  increasing  shock  strength,  the  shock  moves  closer  to  the  nozzle  and 

eventually  vanishes  in  the  nozzle,  i.e.,  x  ->-0;  this  occurs  when  no 

s 

pair  of  Mach  numbers  can  be  found  to  fulfill  the  relation  (26).  In  this 
case,  to  be  discussed  later,  a  shock  may  still  occur  in  the  nozzle. 

For  the  supersonic  branch  M^  =  f(x)  results  from  integrating 
equation  (2U)  with  the  boundary  condition:  M  =  M  at  x  =  0,  where  M^ 

(>  l)  is  given  from  the  area  ratio  of  the  nozzle  by: 
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i«K 

X 
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Figure  3.  Outer  Solution  Variation  of  Mach  Number  in  Gas  Tube 


“  -  f- “  m"  [^tr  (1  +  1rMp2)) 

A  <»  P 


V+l 

2?y<F 


The  subsonic  branch  =  g(x)  is  obtained  by  integrating  equation  (2k) 
with  the  boundary  condition 

x  =  1:  M  =  M  =1  for  p  >  p 
e  *e  rc 

For  pg  =  pc>  neither  Mg  nor  is  known  beforehand,  but  must  be  deter¬ 
mined  by  an  iterative  process.  First,  equation  (2U)  is  integrated  for 
various  values  of  Mg  <  1,  thus  giving  a  set  of  value  pairs  ML^  *=  F(Mg). 
Another  relation  between  and  is  readily  found  from  the  continuity 
equation  and  the  thermal  equation  of  state 


Pn<xa> 


pnus! 


Here,  the  index  II  designates  the  quantity  on  the  subsonic  side  of  the 

shock  at  the  location  x  .  For  the  ”no  shock”  condition  x  =  0  and 

s  s 

P-^O)  =  Pp*  In  any  ca8e’  Pn(xs)  can  ^e  uniquely  related  to  p^  and, 
thus,  to  pg  as  will  be  shown  below.  The  ratio  of  static  temperatures 
follows  from  the  energy  equation: 


e  tot 


1  +  Hi. 

x  2 


0II  tot  1  +  M  2 
2  e 


This  relation  is  valid  for  any  station  x,  i.e.,  for  any  index,  whether 
a  shock  occurs  or  not.  Since  pc  and  pg  are  given  for  each  step  in  the 
parameter  t,  equation  (28)  is  another  relation  of  the  form  =  G(Mg). 

The  two  relations  =  F(Mg)  and  =  G(Me)  can  be  solved 

numerically  for  and  Mg. 

With  Mg  known  the  other  variables  with  the  subscript  e  are  readily 
found.  These  variables  are  used  in  obtaining  the  inhomogeneous  terms 
in  the  differential  equations  for  the  cavity  flow  in  Section  III. A  and 
are  functions  of  the  parameter  t. 
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The  state  with  the  index  e  relates  to  the  state  on  the  subsonic  side 
of  the  shock  by: 


M  { 
e  t 


e 


u  M 
e  e  _e_ 


[w. 


h_  W 

pii  hi  °e 


The  normal  shock  relations  connect  state  I  with  state  II  and  are  given 
here  for  completeness : 

pn  2y  -  (y-1) 
p7"  Y  +  1 

(31) 


PII  = 
Pt 


(y+i)  Mj 
(y-1)  MjS 


The  state  I  is  connected  to  state  p  hy  equations  analogous  to  equation 
(30). 


*I(V 

M 

T> 

pp 

'  “A1 

v*.> 

W 

u 

P 

M 

P 

pI(xs> 

*  ■  ■  ■ —  ■ 

PP 

pp 

P 

qtUT 


©l(xg)/0p  is  given  hy 


31* 


(33) 


0I(xsJ 


0 

P 


0 


I  tot 


0 


L +J? 


If  no  shock  occurs  in  the  duct  the  state  e  or  any  other  state  x  is 
related  to  state  p  immediately  by  equations  of  the  form  (30).  These 
equations  provide  all  the  variables  at  the  entrance  of  the  cavity  for 
all  times  except  for  very  early  times  where  the  quasisteady  approximation 
is  invalid. 


Keeping  in  mind  that  the  very  short  nozzle  is  considered  frictionless, 
the  variables  with  the  index  g  and  p  are  related  by  the  following  if  Mp  >  1: 


0/0  =1+  ^  M  2 
g  P  2  p 


p  /p  =  (1  +  M  2)Y/<Y-1) 
g  P  2  p 


p  /p  = 
g  P 


2  p 


(34) 


=  Pp/P 
=  Mp  (YPp/Pp) 


If  no  shock  can  be  found  in  the  duct  then  Mp<  1,  but  a  shock  may 

still  occur  in  the  nozzle.  This  is  the  case  where  Mp  is  too  large  to  give 

a  purely  subsonic  flow  in  the  nozzle.  The  limiting  may  be  computed 

from  equation  (27)  which  also  allows  a  subsonic  solution.  For  any  Mach  number 

1  >  M  >  M  , .  there  occurs  a  shock  in  the  nozzle.  The  density  p  may 
p  p  lim  p 

be  computed  from 


PP  Pg  (Y+1) 


i/(y-i) 


(ojM  ) 
P 


-1 


Y-l  2  -1/2 

(1  +  ^2~  Mp  )  (35) 
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The  condition  p  >  p  may  occur  and  in  this  case  there  is  back  flow 
c  g 

from  the  cavity  to  the  barrel.  This  back  flow  occurs  with  p  -  p  «  p 

eg  c 

so  that  the  effect  of  pressure  on  the  density  may  be  neglected,  i.e., 
p  j.  fn  (p).  However,  the  effect  of  temperature  on  the  density  may  have  to 
be  considered. 

From  equation  (2i)  and  conservation  of  mass  and  momentum  applied  across 
the  discontinuities,  there  follows  for  the  backward  flow  in  accordance  with 
Figure  4* 


ue" 


2  .  f ,(?c ZlsL 


L  e 


1)  +  to2  [1  +  (—  -  1)  ]  +  — 

°p  5 


(38) 


Here  a  and  a  are  contraction  coefficients.  These  are  a  strong  function 
e  p 

of  the  Mach  number  with  ag  -  0.62  for  incompressible  flow  (M  =  0)  and 

c»e  =  1.0  for  M  =  1.  eip  is  also  a  function  of  m.  For  the  compressible  flow 

the  contraction  coefficients  have  seen  set  equal  to  1.** 

The  coefficient  f  includes  losses  due  to  bends  as  well  as  the  friction 
in  the  duct  and  is  given  by 

f  =  [5/(4S)}  (Kj,  +  Kg)  (39) 

The  pressure  at  station  e"  is  related  to  the  pressure  in  the  cavity  by 

P, 


-  pe»  =  f  V2  [(^“-  1)2  +  X] 

a  w 


(4o) 


with  0^,=  0c,  pe„  =  pc,  and  where  pgfl  pg,  etc.,  in  equations  (l),  (2),  (3). 


The  later  phases  of  the  forward  flow  may  also  be  computed  using 
p  £  fn  (p)  but  p  =  fn  (e). 


*  Derivation  in  Appendix  B 

**For  p  ^  fn  (9),  the  assumption  p  f  fn  (p)  gives  incompressible  flow. 
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The  coefficients  a  ,  ag,  a  ,  and  Kp  and  are  listed  in  Table  3. 


D.  The  Inner  Problem 

As  mentioned  earlier  the  solution  given  in  Section  III.C  is  not 
uniformly  valid  in  t  and  fails  for  small  t.  This  singular  behavior  is 
a  result  of  loss  of  the  time  derivative  in  the  differential  equation  by 
the  expansion  procedure.  The  nonuniformity  occurs  near  the  line  t  -  0 
and  in  order  to  investigate  this  region  a  magnified  variable  T  is 
introduced.  The  variable  x  remains  unchanged  since  the  quasisteady 
solution  behaves  regularly  here.  The  dependent  variables  in  the 
earlier  times  are  of  the  same  order  as  for  the  later  times  and  also  re¬ 
main  unaltered.  Thus, 

T  =  t/6(e)  U  =  u  a  =  p 

X  =  x  P  =  p  H  =  h 

A  =  a 


Introducing  these  variables  into  the  differential  equations  (17),  (l8), 
and  (19)  gives 


e  an 

6(e)  3T 


a(nu) 
a  x 


_L_  l+nl+l  3P 
6(e)  aT  ax  n  ax 


0 


(1*3) 


e  a 
6leT  aT 


e  ap 
6(e)  3T 


0 


where  we  have  dropped  the  friction  and  heat  transfer  terms  in  accordance 
with  the  discussion  in  Section  II.  In  the  limit  e  +  0,  e/5(e)  0 
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results  again  in  the  loss  of  the  highest  derivative,  while  e/5(e)  » 
when  e  •*  0  leads  to  a  solution  which  cannot  be  matched  with  the  outer 
solution .  Thus ,  we  must  have 


-*■  const,  as  e  0 


Without  loss  of  generality  the  constant  is  taken  equal  to  thus,  we 
have  6(e)  =  e,  and  we  get  the  differential  equations 


an  a(au)  = 

sf  +  sr  a 


0 


1 

n 


3P 

3X 


=  0 


3_ 

3T 


+ 


U 


3_ 

3X 


(H 


with  the  initial  conditions 


T  «  0,  X  >  0:  P  =  Pi; 


U  =  0 


(44) 


X  <  0:  P  =  P  ;  A  =  A  ;  U  =  0 

g  g 

The  boundary  conditions  for  T  >  0  are  those  existing  at  the  edge 
of  the  region  of  nonuniforaity . 

X  =  0,  T  >  0:  P(0,T)  =  p(0,t),  A(0,T)  =  a(0,t) , 

t  -*■  0  t  -*  0 

U(0,T)  =  u(0,t) 
t  0 

For  X  ~  1  no  boundary  condition  is  applied  as  long  as  Ug  >  Ag;  and  for 

U  <  A  the  condition  is 
e  e 

T  >  0,  X  =  1:  P(1,T)  =  pc(t)  =  pt 

t  ■+■  0 
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There  is  no  continuous  solution  satisfy  * ng  these  initial  and  boundary 
conditions  and  the  actual  solution  contains  surfaces  of  discontinuities. 
A  solution  is  readily  constructed  using  shock  wave,  contact  front, 
uniform  region  of  flow,  unsteady  expansion  and  steady  expansion.  The 

problem  is  that  of  a  shock  tube  having  an  area  change  at  the  diaphragm 

.  ,,  8,9,10 

station.  *  * 


Figure  5  shows  an  X,T  diagram  of  the  flow  in  the  customary  shock 
tube  notation.  Since  one  is  interested  in  the  influx  of  powder  gas 
into  the  cavity  one  may  neglect  the  gas  between  shock  and  contact  front, 
which  i?,  of  course,  the  gas  (air)  contained  in  the  tube  at  time  t  <  0, 
and  is  not  powder  gas. 

With  U 0  =  U,  across  the  contact  front  we  set  U  =  0  for  T  <  U 

In  the  region  3  we  have  constant  flow  properties.  Thus, 

Ue  =  U3  =  U2  =  const.  for  l^-1  <T<  [Ug^-  (y-l)  R]"1 


where  2R  =  Up  +  A^  =  Const.  (=  U^,  +  A^)  is  the 

Riemann  invariant.  The  pressure  across  the  contact  front  is  equal  so 
that  P2  =  P  . 

P  follows  from  the  shock  relation  (31).  For  M  »  1  we  have 
e  s 


Pe/Pi  =  P./P,  •  P-  M 
Y+l 


<Ms= 


and  rewriting  this  expression  in  terms  of  Ugi 


P2/P 


1 


Y  (Y  +  1) 
2 


(U5) 


8 


Ug  may  be  found  by  numerically  solving 


Pjo) 


.  ...  u_  2 

JL_  *  (_2) 

P  lA  ' 

1  1 


/.  A 


v  (v  +  l) 

2g 


l  - 


n  _  III 

2  y-1  2y 

r To)  2  g 

g 


_  2l_ 
Y-1 


A 

Here,  y  -  1.4  and  the  gain  factor  g  is 


S=< 


r 

ft>  “ 

,  .  2  . 

1/2 

'2  +  (y-1)  Mr]| 

U  +  (y-1)  m3,s. 

2  | 

(46) 


where  M.^  =  M  from  equation  (27).  Flow  variables  at  state  3'  are  found 

^  Jr 

from  equation  (34). 


Vap  -  w 


V  (y+1) 
2 


/M2  P,  1/Y 


(hT) 


It  is  emphasized  that  constant  flow  properties  in  this  region  are  a 
result  of  neglecting  heat  transfer  and  friction  and  also  shock 
reflections;  actually,  these  variables  depend  on  Reynolds  number,  geometry 
and  time. 

The  flow  properties  in  the  expansion  fan  depend  on  T  for  X  =  1. 

Thus,  for: 

[u2  -  (y-1)  Rj  <  T  <  [V  (y-1)  r| 


we  have 


U(1,T)  =  $  2R  .  ^  i 


P(1,T)  =  P, 


f-i-  2=i  (2R-i-)l  Y_1 

[A  .  y+1  '  n  t'J 


(48) 


a(i,t)=  iilfa,  -i) 

•y  +  1  '  T  / 
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For  time 


Jiy  (y-D  b]"1  <  T  < 


the  flow  properties  are  constant 

U(1,T)  =  Up  =  U3.  A(1,T)  =  Ap  =  hy 

(49) 

P(1,T)  =  Pp  -  Vy 

E.  The  Composite  Solution  for  the  Flow  in  the  Duct 

The  analysis  described  so  far  has  yielded  the  first  terms  of  two 
complementary  asymptotic  expansions.  In  general  the  two  asymptotic 
expansions  valid  for  small  and  large  times  respectively  have  a  domain 
of  overlap,  in  which  both  are  valid.  This  is  the  case  whenever  Van 
Dyke,s"asymptotic  matching  principle"  is  valid.  This  principle  applied 
here  amounts  to  setting  the  first  term  inner  expansion  of  the  first 
term  outer  expansion  (the  outer  expansion  rewritten  in  inner  variables 
and  expanded  for  small  e  truncated  to  the  first  term)  equal  to  the 
first  term  outer  expansion  of  the  first  term  inner  expansion.  This 
principle  leads,  for  the  velocity ,  to  the  condition 

U(X,«)  =  u(x,0)  (50) 

with  corresponding  expressions  for  the  other  dependent  variables.  We 
note  that  from  equation  (49) 

U(l,»)  =  Ug(»)  =  by  =  Up 

For  inviscid  and  nonheat-conducting  flow  the  solution  to  equation 
(21)  is  trivial  and  gives 


u(x:  t).  .  =  u  (t).  .  =  u  (t).  .  =  const. 

1  mvisc  p  mvisc  e  mvisc 

where  the  constant  is  a  function  of  the  parameter  t. 


For  t  =  0  ve  have 


u(*-0)invisc 


u 


p 


invisc 


(0) 


suM 


so  that  the  matching  condition  is  fulfilled  for  inviscid  flow.  In  this 
case  the  two  solutions  may  be  combined  to  yield  one  composite  solution, 
which  then  is  uniformly  valid  over  the  whole  range  of  t. 


Using  the  multiplicative  composition  of  reference  11,  one  finds  the 
composite  expansion  as  the  outer  expansion  multiplied  by  a  correction 
factor  consisting  of  the  ratio  of  the  inner  expansion  to  its  outer 
expansion. 


u 


corap. 


(51) 


Since  in  the  present  case  the  outer  solution  has  been  computed  on 
the  basis  of  viscous  heat-conducting  flow,  matching  of  the  inner  and 
outer  solutions  is  not  immediately  possible  for  x  >  0.  For  x  =  0  and 
supersonic  flow  at  the  entrance  of  the  duct,  matching  of  the  inner  and 
outer  solution  according  to  equation  (50)  is  effected;  for  subsonic  flow 
matching  is  not  possible,  reflecting  the  fact  that  the  disturbances  due 
to  friction  now  affect  the  entrance  condition  in  the  duct.  Failure  to 
match  according  to  (50)  of  course  simp]y  represents  the  fact  that  the 
inner  and  outer  solutions  axe  different  approximations  in  the  parameter  1/Re. 


Formally  the  composite  solution  of  the  form  ( 51 )  may  be  used  here 
also.  The  effect  of  friction  and  heat  transfer  on  the  velocity  may  be 
assumed  accounted  for  in  the  following  way: 

U(l,T)visc  =  B  U(1,T)  (52) 

where  the  factor  8  will  depend  on  Reynolds  number,  geometry,  Mach  Number 
and  pressure  ratio  and  on  time.  The  composite  solution  will  be 


«  =  ud.t)  SW 

comp. ,  e  U(l,«)  8(“) 


(53) 


For  t/e  <  (ly)’1  ,  ucompf 


=  0  because  U  =0.  This  time 


t  =  e  x  (UgB)  =  i  (UgB)  is  of  course  the  delay  time  between  rise  of 


pressure  at  the  port  and  rise  of  pressure  in  the  cavity.  This  time  has 
been  observed  experimentally ;  it  determines  3,  the  average  value  of  B 
in  the  range  0  <  t  <  £  (UgB)  ^  .  The  experimental  values  of  B  scatter 

between  1  and  0.5. 

For  t  =  “  one  can  find  from  the  matching  condition,  when  applied  to 
viscous  flow,  B  {“)  as 


B  (»)  =  u(l,0)  /  U(l,«), 

and  this  leads  to  a  numerical  value  of  B  ~  0.35.  Since  B(t)  is 
expected  to  vary  quite  rapidly  near  t  =  0,  it  is  assumed  that  at  the  end 
of  the  first  interval,  i.e.,  at  t  =  I  (U^B)  \  B  has  essentially  reached 
the  value  for  t  -»•  «;  thus  we  set  the  ratio  B  (t/e)  /  3  (»)  -  1. 

The  ratio  3  (t/e)  f  B  (»)  can  also  be  considered  an  adjustable  pa¬ 
rameter  with  which  the  computation  could  be  adjusted  so  as  to  give  best 
agreement  with  experiment.  This  is  not  done  here  since  this  procedure 
could  disguise  other  phenomena.  It  is  not  possible  to  correctly  account 
for  the  friction  and  heat  transfer  at  the  very  early  times,  as  has  been 
pointed  out  before,  and  some  sort  of  plausible  assumption  has  to  be 
made.  The  basic  difficulty  is  that  the  equations  (13)  and  (l4)  with 
the  quasisteady  friction  and  heat  transfer  coefficients  are  not  appli¬ 
cable  in  this  early  phase;  and  the  treatment  of  this  problem,  say,  by 
the  method  of  characteristics  would  of  course  not  remove  this  difficulty. 

With  the  assumption  that  the  viscous  correction  factor  for  velocity 
B(t/e)  /  B(°°)  -  1  for  t  >  £  (U^B)  ^  and  that  the  corresponding 
correction  factors  for  the  other  variables  are  also  approximately  unity , 
che  effect  of  friction  appears  only  in  the  delay  time  t  =  £  (U^B)  , 

where  it  can  be  observed  experimentally,  and  in  the  other  interval 
boundaries,  where  it  cannot  be  directly  observed. 

The  composite  solution  in  the  appropriate  intervals  is  then 


u  =0, 

conip ,  e 

for  Q  <  t  <  StUgB)"1 


u  =  u(£;  t)  0  /U  , 

uomp,e  2  3 

for  l  (U2B)  “1  <  t  <  £  /  [U2  (*—)  -  (y-1)  R ]0 


u  =  u  Hit)  (*=f  2R  +  -)  U„._1 

comp,e  Y+l  Y+l  £  3 


for  4  /  [Ug  -  (y-1)R  ]  6  <  t  <  £  /  [0  -  (y-1)R]  B 

u  =  u  (£;t) 

comp,  e 

for  l  /[U  '  ,  (*~4  -  (y-1)  R]  B  <  t  <  »  . 

It  has  "been  assumed  here  for  simplicity  that  B  for  the  Riemann  in¬ 
variant  is  the  same  as  6  for  the  velocity.  This  assumption  is  quite 
acceptable  within  the  assumptions  made  already  and  can  only  result  in  a 
slight  shift  in  the  interval,  boundaries.  It  should  be  remembered  that 
the  time  after  which  the  effect  of  the  first  phase  has  completely  died  out 
is  t  =  £  /  [U3<.(^~)“  (y-1)  R]  8,  and  that  this  time  is  still  small 
compared  to  r. 

The  composite  solution  for  the  pressure  in  the  corresponding  intervals 


p  =  p  (£;t)  P, /P,. 

comp,  e  e  13 


p„„  „  =  p  (£;t)  p9/p,. 

comp ,  e  e  2  5 


p  =  p  {1 

*comp ,  e  e 


;t)  [—  ~r-  (2  R  -  — ) 

.V  t  J 


2y/(y-1) 


p  =  p  (£;t) 

rcomp ,e  Q 


Similarly  for  the  sound  speed 


a 

comp*  ,  e 


l&  (£;t)  (P2/P3J(y"1)/2y 


(56) 


The  composite  solutions  given  by  equations  (51*),  (55),  and  (56) 
provide  ug,  pe,  and  ag  for  equations  (l),  (2),  and  (3).  The  other 
dependent  variables  at  the  cavity  entrance  (location  e)  may  be  computed 
from  these  quantities. 

IV.  COMPARISON  WITH  EXPERIMENT 

A  comparison  between  theory  and  experiment  was  made  on  the  M-l6 
Rifle.  Input  data  for  the  theoretical  computation  consist  of  the 
measured  time  variable  pressure  in  the  barrel  at  the  port  location,  the 
physical  properties  of  the  powder  gas,  and  the  dimensions  of  the  gas 
system.  The  physical  properties  of  the  powder  gas  are  listed  in  Table 
1.  They  represent  estimates  based  on  the  interior  ballistics* 


Table  1.  Properties  of  Powder  Gas 


=  1.7**  x  103  (m/sec)2/°K 
=  1.38  x  103  (m/sec)2/°K 


=  0.831  x  10-1  (kg  m)/(sec3  °K) 
=  0.U00  x  103  (ra/sec)2/°K 


=  1.26 
=  1.2U 

=  U.80  x  10  kg/(m  sec) 


■^Supplied  by  Mr.  R.  Geene,  of  Interior  Ballistics  Laboratory 


computation  of  the  composition  of  the  powder  gas.  The  dimensions  of 
the  gas  system  are  summarized  in  Table  2.* 


A 


min 


Table  2.  Dimensions  of  Gas  System 

.  -k  2 
=  1.269  x  10  m 

=  0.6605  x  10"5  m2 

=  0.1*383  x  10“5  m2 

=  0.6605  x  10"5  m2 

=  0.6533  x  lO"5  m2 

=  0.2900  x  lO-2 

=  0.3600  m 

=  O.U366  kg 


m 


-2 


m 


=  0.11*1*2  x  10' 

=  0.7600  x  10"  nT 

=  0.7600  x  l<f2  m 
-2 

m 


D 
i 

% 

r 

v 

V  . 
ci 

*BE 

=  0.6158  x  10 
(e/D)  «  0.015 

Since  it  is  necessary  to  specify  another  variable  of  state  at  the 
port,  the  temperature  at  the  port  was  computed  according  to: 

0(t)/0(O)  =[p(t)/p(0)]  <?**!}/?  C57} 

Here  y  is  not  the  ratio  of  specific  heats,  but  an  exponent  determined  by 
fitting  the  temperature  variation,  known  from  the  standard  interior 
ballistic  treatment  of  the  weapon,  to  the  pressure  variation  according 
to  equation  (57).  The  standard  interior  ballistics  treatment  breaks 
down  after  the  bullet  has  left  the  barrel;  the  temperature  distribution, 
however,  is  needed  for  the  computation  after  this  time,  and  is  extra¬ 
polated  according  to  equation  (57). 


♦Supplied  by  Mr.  M.  Werner  of  Interior  Ballistics  Laboratory. 
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Probably  the  largest  ■uncertainties  lie  in  the  various  los3  factors. 

For  the  friction  factor  the  value  for  clean  commercial  steel  pipe  has 
12 

been  used.  Compressibility  effects  on  the  friction  are  apparently 

7 

quite  small,  at  least  for  hydraulically  smooth  pipes  -where  Frossel  did 

not  find  any  dependence  on  the  Mach  number  of  the  flow.  For  rough  pipes 

such  as  the  commercially  clean  pipe  used  as  duct  in  the  M-l6  rifle, 

(relative  roughness  i/D  -  0.015)  the  dependence  of  the  friction  factor 

on  Mach  number  is  probably  more  pronounced,  and  will  result  in  a  slightly 

smaller  friction  factor.  Since  during  the  time  t  the  flow  in  the  duct 

ranges  from  supersonic  to  subsonic  and  incompressible  flow,  the  friction 

factor  could  be  introduced  as  depending  on  the  Mach  number.  This  has, 

however,  not  been  done;  the  friction  factor  is  assumed  to  be  constant 

at  the  value  for  incompressible  flow.  It  may  be  mentioned  that  the 
2 

experimental  work  was  conducted  with  a  used  rifle  showing  evidence  of 

deposits  in  the  interior  of  the  duct.  X-ray  pictures  of  ducts  in  the 

M-16  rifle  show  indeed  substantial  amounts  of  deposits  depending  on 

13 

powder  type  and  usage  of  the  rifle.  The  deposits  are  so  located  as  to 
cause  little  additional  resistance  in  forward  flow  where  stagnation 
pressure  is  high,  but  in  the  backward  flow  they  may  cause  a  considerable 
increase  in  losses. 

The  losses  due  to  bends  are  taken  into  account  by  loss  coefficients 
as  is  done  in  incompressible  flows.  The  loss  coefficients  are  taken 
from  reference  12.  The  losses  apply  to  miter  bends,  which  are  most 
representative  of  the  bends  in  the  gas  duct  system.  Again,  compress¬ 
ibility  effects  probably  tend  to  decrease  the  losses  slightly. 

The  area  reduction  at  the  port  has  been  treated  as  a  frictionless 
nozzle  for  supercritical  pressure  ratios  and  as  an  abrupt  area  change 
for  subcritical  pressure  ratios  as  is  explained  in  Section  III.  In  the 
Addendum  the  analysis  has  been  extended  to  include  treatment  of  the  area 
reduction  as  an  abrupt  area  change  e^en  for  the  case  of  supercritical 
pressure  ratios.  In  this  case  supersonic  flow  cannot  be  reached  in  the 
duct.  It  is  not  possible  to  decide  which  of  the  two  types  of  flows  is 
actually  established,  because,  as  will  be  shown,  the  effect  is  easily 
obscured  by  changes  in  the  friction  factor. 
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For  incompressible  flow  the  area  reduction  as  well  as  the  sharp 

edge  entrance  causes  a  decrease  in  mass  flow  which  is  accounted  for  by 

12  ll* 

contraction  coefficients.  * 

Table  3  gives  a  summary  of  all  Lhe  loss  coefficients  usee  in  the 
computations. 

Table  3.  Loss  and  Flow  Coefficients 

a  =  0.62  Kg  =  5.0U 

a?  «  0.6U  Kg  =  6.95 

a  -  0.62  f  =  0.02lll 

g 

2 

Figure  6  shows  a  comparison  between  computed  and  experimental 
pressure  variations  with  time  in  the  bolt  cavity.  It  is  seen  that 
agreement  between  experiment  and  theory  is  quite  satisfactory.  The 
theoretical  pressure  is  rising  somewhat  more  slowly  than  the  experi¬ 
mental  value.  This  could  be  due  to  a  friction  coefficient  somewhat 
too  large  in  the  early  compressible  phases  of  the  flow.  However,  the 
fast  experimental  rise  could  more  plausibly  be  explained  by  the  fact 
that  there  is  a  component  of  dynamic  pressure,  resulting  from  the  gas 
motion  in  the  barrel  itself,  acting  on  the  port,  thus  effectively  in¬ 
creasing  p  ;  this  would  result  in  a  larger  influx  in  the  bolt  cavity, 
g 

giving  a  correspondingly  faster  pressure  rise.  Indeed,  this  is  sug¬ 
gested  by  the  fact  that  barrel  material  is  "washed  out"  downstream  of 
the  bore  hole  comprising  the  port.  This  same  mechanism  would  also  re¬ 
duce  the  backward  flow  since  the  effective  p  is  larger  than  indicated 

g 

by  the  static  pressure  measurement.  This  is  also  borne  out  by  the 
comparison.  It  is  seen  that  the  predicted  pressure  decreases  faster 
than  the  experimental  pressure;  this  is  caused  by  too  large  a  back  flow. 
(See  marker  on  Figure  6.)  As  has  been  pointed  out  above,  the  resistance 
due  to  deposits  is  also  expected  to  be  larger  in  the  back  flow  than  in 
the  forward  flow,  and  this  has  not  been  accounted  for  in  the  model. 


Figure  6.  Comparison  of  Experimental  and  Theoretical  Bolt  Carrier 
Pressures  for  M-16  Rifle  (Round  5?) — Nozzle  at  Port 
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Th_  measurements  were  made  on  a  recoiling  gun.  This  means  that  the 
piston  displacement  x^  is  measured  from  an  accelerated  reference  system. 
As  a  result  there  is  an  apparent  force  on  the  bolt  carrier  which  opposes 
the  bolt  carrier  motion.  The  measurement  of  gun  motion  indicates  that 
the  gun  acceleration  has  dropped  to  nearly  zero  when  the  pressure  starts 
to  rise  in  the  bolt  cavity,  so  that  only  the  early  part  of  pressure 
history  could  be  affected.  It  is  worth  noting  that  the  net  effect  of 
the  acceleration  would  result  in  a  faster  pressure  rise  than  without  the 
acceleration.  The  fast  experimental  pressure  rise  could  indicate  that 
the  gun  is  still  accelerating  when  the  pressure  starts  rising  in  the 
cavity. 

In  this  connection  it  should  be  mentioned  that  in  the  numerical 
computation  the  friction  forces  ®  on  the  bolt  carrier  were  assumed  zero. 
These  are  indeed  expected  to  be  small  compared  with  the  inertia  force. 

The  venting  area  has  been  assumed  to  be  the  nominal  area  of  the 
venting  holes.  Inspection  of  a  number  of  guns  indicates  that  this  area 
may  not  be  "cleared"  completely  by  the  bolt  due  to  manufacturing  toler¬ 
ances.  This  would  affect  the  tail  of  the  pressure  curve  only  and  will 
not  affect  the  motion  of  the  bolt  carrier. 

Figure  7  shows  the  comparison  of  the  experimental  pressure  trace 
with  the  computation  where  the  area  change  is  treated  as  a  discontinuity. 
The  computation  based  on  the  friction  factor  f  =  0.02U1  (same  as  in 
Figure  6)  gives  a  somewhat  lower  pressure,  but  decreasing  this  friction 
factor  by  25#  gives  essentially  the  same  curve  as  shown  on  Figure  6. 

Since  the  friction  factor  is  not  believed  to  be  known  to  better  than 
25$,  it  is  not  really  possible  to  make  a  choice  between  the  two  models. 

V.  CONCLUSION 

The  analysis  presented  here  gives  the  first  term  of  an  expansion 
in  the  parameter  e,  which  is  the  ratio  of  acoustical  transit  time  to 
the  characteristic  time  of  pressure  variation  impressed  on  the  flow  in 
the  actuation  mechanism  of  the  automatic  weapons.  Automatic  weapons,  by 
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Figure  7.  Comparison  of  Experimental  and  Theoretical  Bolt  Carrier 
Pressures  for  M-l6  Rifle  (Round  52)--Area 
Discontinuity  at  Port 
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necessity,  are  designed  so  that  this  parameter  is  small  in  order  to 
assure  quick  responsive  operation  of  the  extracting  and  loading  mechanism,, 
Thus,  although  this  analysis  has  only  teen  tested  for  the  M-l 6  rifle,  it 
is  expected  to  apply  to  other  weapons  and  weapon  systems  also. 

Most  of  the  automatic  weapons  actually  have  considerably  smaller 
values  of  e  because  the  length  of  the  duct  is  very  much  shorter  than  it 
is  the  M-l6  rifle  and  many  weapons  have  no  duct  at  all.  '  For  these 
weapons  the  present  analysis  is  also  applicable,  if  the  inner  solution 
is  completely  discarded  and  the  flow  is  treated  as  frictionless. 

An  approximate  treatment  for  the  M-l6  rifle  is  obtained  when  the 
volume  of  the  duct  is  added  to  initial  volume  and  the  inner  solution  is 
discarded.  The  heat  losses  due  to  heat  transfer  in  the  pipe  have  to  be 
considered  and  may  be  accounted  for  globally  in  equation  (3).  This  is 
the  approach  used  in  reference  2.  In  this  treatment  the  initial 
condition  is  not  fulfilled  and  the  solution  is  not  valid  for  small  times; 
however,  bolt  velocity  and  bolt  trajectory,  which  depend  on  integrals 
of  the  pressure  distributions,  are  well  predicted. 

Most  interesting  would  be  the  extension  of  this  analysis  to  include 
second  order  terms.  The  characteristic  time  r  does  not  appear  explicitly 
in  the  first  order  solution;  however,  it  should  not  be  concluded  that  t 
is  an  "artificial"  time.  Actually,  z  will  explicitly  occur  in  the  second 
approximation. 
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ADDENDUM 


TREATMENT  OF  PORT  AS  AREA  DISCONTINUITY 


The  assumption  used  in  the  body  of  the  report  treating  the  area 
constriction  at  the  port  as  a  frictionless  supersonic  nozzle  is  an 
artifice,  especially  since  conditions  are  encountered  which  require  a 
standing  shock  in  the  nozzle. 

It  may  "be  more  appropriate  to  treat  the  area  constriction  as  an 
area  discontinuity.  In  this  case  supersonic  flow  cannot  he  reached  in 
the  duct.  Instead  is  known  from  the  integration  of  equation  (24). 
Then  the  quantities  at  the  throat  (*  quantities)  may  he  computed  from 

equations  (36)  and  (37).  Equation  (36)  is  only  applicable  if  p  <  p 

*  *  crlt 
while  for  p  >  p  ^  we  have  M  =  1.  It  may  he  shown  that  the  latter 

condition  is  equivalent  to  the  condition 


“  Y  + 


where  G  is  given  as  G(M  )  with  equation  (36).  This  inequality  is 

P 

always  fulfilled  for  =  1.  Consideration  of  momentum 

p  A  [1  +  y  M  2]  +  p  (A  -  A)  =  p  A  [1  +  v  M  2] 

P  P  P  P 

* 

then  shows  that  o  -*•  0  for  M  =  M  =  1;  i.e.,  the  pressure  on  the  surface 
#  P 

Ap  -  A  is  zero.  (The  control  volume  for  application  of  the  momentum 
equation  is  shown  in  Figure  8.) 

The  inner  solution  also  has  to  he  modified  to  account  for  the 
discontinuous  area  change.  If  the  flow  between  throat  and  station  p 
is  considered  steady,  then  M  =  1  and  the  shock  velocity  follows  from 
solving  numerically 

_  £l_ 

1/2  m2  y-1 

hw/*tmu [Jr*?-  t 


CONTROL  VOLUME 


Figure  8.  Sketch  for  Compressible  Flow  Treatment 
of  Area  Discontinuity  at  Port 
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The  inner  solution  is  shown  in  the  X,T  diagram  of  Figure  9*  The 
relations  (1*5),  (U7)  and  ( U8 )  are  valid  here  also.  The  conditions  at  p 
now  follow  from  (37)  consistent  with  the  assumption  of  steady  flow 
"between  throat  and  station  p.  It  is  noted  that  the  expansion  fan  now 
"reaches  hack"  to  station  p  and  the  inner  solution  affects  the  flow  in 

the  duct  even  for  T  -*■  ». 

« 

The  condition  M  =  M  =1  also  makes  the  pressure  on  the  area 

«  P 

A  -  A  equal  to  zero,  as  has  been  noted  before.  It  is  worth  pointing 
out  that  this  steady  flow  condition  is  reached  as  the  limit  of  the 
unsteady  flow.  The  unsteady  flow  would  result  from  the  reflection  of 
the  expansion  fan  from  the  surface  Ap  -  Jt;  and  for  A  4  Ap  the 
reflection  will  be  similar  to  the  reflection  of  an  expansion  from  a 
solid  wall.  It  is  well  known  that  vacuum  may  result  on  the  wall  for  a 
sufficiently  large  pressure  ratio  through  the  expansion  fan. 
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*  p 


Figure  9.  X,  T  Diagram  for  Inner  Solution — Area 
Discontinuity  at  Port 
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APPENDIX  A 


DERIVATION  OF  EQUATION  (36) 


From  the  conservation  of  mass 


*  *  *  ,*  ,  ,1/2 

p  a  m  =  p  a  m  (e/e  )  ' 

p  p  p  p 


and  momentum 


#  *  ,  *2,  #  *  ,  .  2, 

pA(l+yM)+pA  (to-1)  =  p  A  (l  +  y  M  ) 


there  follows 


(M  /M)  (0/6  i1^2  (to  +  y  M2)  =  1  +  y  M  2 
P  P  P 

*  . 

Replacing  0/0  by  the  energy  equation  (37)  and  solving  a  resulting 

p  *2 

quadratic  equation  for  M  lead  to  equation  (36). 


APPENDIX  B 


DERIVATION  OP  APPROXIMATE  FLOW  RELATIONS 


The  relation  (38)  for  the  flow  velocity  u  at  station  e"  (Figure  , 
•where  the  fluid  again  fills  the  total  cross-sectional  area  of  the  duct, 
may  he  derived  hy  first  considering  that  the  flow  ohey3  Bernouilli's 
equation  from  station  c  to  station  e  * .  Thus 

Pc  ”  Pe.  =  Pue«/2  (B-1) 

3etween  station  e*  and  e"  the  spreading  out  of  the  jet  occurs  under 
losses,  and  application  of  conservation  laws  for  mass  and  momentum 
leads  to 


pe'  “  peu  “  pe  V  [l  "  (l/ce)] 


(B-2) 


where  i3  the  contraction  coefficient  at  station  e.  Combining  (B-l) 
and  (3-2)  and  applying  conservation  of  mass  gives; 

2 


pe  V 

Pc  “  Pe«  =  — ” 


a 


il2  +i 


(B-3) 


Between  stations  e"  and  p  integration  of  the  momentum  equation  (21) 
yields 

.1 


p  u  (u  —  u  „)  +  p  -  p  „  + 
P  P  P  e  p  e 


2  2  .  r 

)  u  f  /  . 

-E — E —  /  I 

2D  J  p 
J  0  • 


dx  =  0 


(B-10 


Combining  (B-3)  and  (B— 1+ ) ,  noting  that  =  pg  =  pc  =  const.,  leads  to 


■d  -  p  = 
c  p 


P„  UQ" 
c  e 


2  r 


[i-  _  1 

l2  .  ,  .  It  f  l 

Ct 

+  1  +  ~E 

l  e 

m 

(B-5) 


The  pressure  difference  between  stations  p  and  g  follows  again  by 
applying  Bernouilli's  equation  between  stations  p  and  p*  and  the  con- 
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servation  laws  between  p'  and  g.  Thus: 


p  u 

■n  —  -n  —  &  ■  ■P 

Pg  pp  2 


/.  » 

/  1  ,  \  2  ' 

\ 

3 

1 

H 

rir-1)  J 

) 

(B -6). 


Combining  equations  (B— 5 )  and  (B-o),  noting  that  u^  =  because  = 
Pe,  gives  equation  (38). 

Equation  (Ul)  can  be  derived  in  an  analogous  manner.  First,  the 
pressure  difference  between  stations  g  and  g,!  (Figure  if)  is 


p  U  it 

p  _pt(=  liiLJL_ 

pg  2 


+  1 


(B-7) 


The  pressure  difference  between  stations  g"  and  p  is 

-I  2  ,  -1  , . 

Pg«  “  Pp  =  »  Pg  «gti  (»  -  1) 

so  that  the  pressure  difference  between  stations  g  and  p  is 

2  r  2 

j~ - 4  +  (to  -  l)2  +  1 


(B-8) 


P  u 

p  _p  =  J5-JL 

g  "P  2 


g 


(B-9) 


For  the  pressure  difference  between  stations  p  and  e 

^  2  2  ? 

Pp  -  ■  Pe  \  [1  -  Vpp,!  +  —7^—  f  J  ^ 


From  (B~9)  and  (B-10)  one  obtains  equation  ( Ul) ,  noting  that  for 
P  i  fn  (p) 

1/p  =  8,  ./( p  9  ). 
tot  g  g 
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